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The dispersion relations for surface plasmon-polaritons propagating in the Voigt geometry in a 
metal-insulator-metal waveguide with a magneto-optically active dielectric medium are derived. The 
symmetry between the upper and lower interfaces is broken by the introduction of the magnetic field; 
the balance between the field distributions on the two interfaces can be controlled by the applied field. 
This control is illustrated by finite-element method numerical simulations of the field distributions 
around a point dipole placed in the centre of the short waveguide; it is shown that both the total 
emission of radiation from the cavity and the distribution of the far-field radiation can be strongly 
modified by tuning the magnetisation of the waveguide. This raises the novel possibility of using 
magnetic fields to control light propagation in nanostructures. 
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PACS numbers: 42.79.Gn,85.70.Sq,73.20.Mf 
I. INTRODUCTION 

Surface plasmon-polaritons (SPPs)^ are electromag- 
netic waves that propagate along metal-dielectric inter- 
faces and can be guided by metallic nanostructures be- 
yond the diffraction limit^'^. They are attractive because 
of the very small length scales over which it is possible to 
localise the electromagnetic fields, giving very wide scope 
for the manipulation of those fields on the nanoscale. 
They usually extend few hundreds of nanometers in the 
surrounding dielectric. They can be further confined and 
guided in metal-insulator-metal (MIM) structures, which 
are building blocks for resonant guide wave networks^ 
and can be used for future on-chip nanocircuits^'^. 

Volume modes in waveguides propagate via multiple 
reflections from the metal surfaces and the energy is 
localised in the core of the structure. Surface plas- 
mon modes, on the other hand, are localised at the 
interfaces; for a structure with two symmetric inter- 
faces they couple to produce symmetric or antisymmet- 
ric mode proflles. The case of metal-insulator-metal 
(MIM) plasmonic waveguides with a non- magnetic dielec- 
tric has been studied previously both theoretically^^^ 
and experimentally In these studies, MIM plas- 
mons were excited by an external light source using, for 
example, in-coupling through slits in one of the metal 
cladding layers, or via the electron beam of a scanning 
electron microscope. 

Placing an emitter inside a MIM cavity is another 
possible route to excite cavity plasmons. Dipolar emit- 
ters placed in MIM slab and slot waveguide structures 
were found to couple strongly to the plasmon waveg- 
uide modes^^. Such a structure can serve as an opti- 
cal nanoantenna^^. Active control of the direction of 
emission of an optical nanoantenna has been achieved by 
electrically controlling its load impedance^^. The electric 
fleld has been shown^^ to modulate the dispersion rela- 
tion of the plasmons propagating at the interface of an 



electro-optically active dielectric and metal. 

A quite different route for active control of plasmons 
is by using a magnetic fleld. The combination of plas- 
monics with magneto-optical materials is particularly in- 
teresting because it introduces a nanoscale interaction 
between light flelds and magnetisation, hence opening 
up the possibility of using either one of these flelds to 
control the other. For example, external magnetic flelds 
could then be used to control plasmonic devices through 
nanoscale analogies of the Kerr and Faraday effects^^; 
alternatively, the light fleld could be used to interact 
with a nanoscale magnetic system. SPPs propagating 
at a single metal-dielectric interface where one or both 
media are magnetic have been investigated theoretically 
and experimentally One notable effect occurs in 
the Voigt geometry (magnetisation perpendicular to the 
propagation direction), where the wave vectors for left 
and right propagation become unequal. There has been 
growing interest in these ideas in recent years, and also 
in experimentally realising active magnetic control of the 
plasmons. For example, in^^ it was demonstrated that a 
d.c, magnetic fleld can be used to modulate the plas- 
mons, and it has also been shown that introducing a 
periodic nanostructure results in an enhanced magneto- 
optical Kerr effect (MOKE) signal^^'^^. For a review on 
magneto-plasmonics see^^. Magneto-plasmonic waveg- 
uides consisting of insulator-metal-insulator (IMI) geom- 
etry with a ferromagnetic metal surrounded by nonmag- 
netic dielectrics^^ or a nonmagnetic metal bounded by 
ferromagnetic dielectrics^^ have been studied for appli- 
cation as active devices in SPP-based optics. Plasmons 
at the interface between a magnetic metal and a photon- 
ics crystal have also been considered^^. 

In this article we consider the case of a MIM waveguide 
containing a magnetic insulator. We study the surface 
plasmon modes of the waveguide, a preliminary inves- 
tigation of which was published in^^. We then employ 
numerical simulation to show the fleld and energy dis- 
tribution inside the waveguide. We excite the modes by 
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placing a dipole in the middle of the structure and show 
that an external magnetic field can switch on and off the 
coupling of the dipole to the surface plasmons. We study 
the out-coupled radiation, the intensity and direction of 
which can also be controlled through external magnetic 
field. 



II. MAGNETO-PLASMON DISPERSION 
RELATIONS 

A. Single magnetic interface 

To understand the behaviour of the SPP guiding modes 
in a cavity we first review the case of a single interface 
separating a metal and a magnetic dielectric, which we 
take to be isotropic apart from the magnetic effects. We 
take the boundary at the plane = 0, with the wave 
propagating in the x-direction and the static magnetisa- 
tion being in the z direction (i.e. in the Voigt geometry). 
For the Voigt geometry only TM modes can propagate on 
this interface^^, so we also take the oscillating H-field in 
the z-direction. The permittivity tensor for a gyrotropic 
medium when M = (0,0, M^) is 



with 












er 
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(1) 



with exy = gMz. Let the half-space ^ > contain the 
magnetic dielectric and ^ < the non-magnetic metal. 
We seek to find waves damped as ^ ±oo in the form 

i^yd = Vkl - w'^ed/c'^, y > 



in the dielectric, and 



(2) 



in the metal, where ed = er -\- ^ly/^r and is the per- 
mittivity of the metal. 

The dispersion relation can be found from the bound- 
ary conditions on the tangential E and normal D fields 
and is given by 

^m^r^ym ^^m^xy^x (^r ~^ ^xy)^yd ~ 0? (3) 

which for the nonmagnetic case exy = results 
in the well-known plasmon dispersion relation kx = 
'^^^ yTe '^-^e ) ' where ko = ou/c is the magnitude of 
the free-space wave vector corresponding to angular fre- 
quency cj, in the region where e^e^ < and + < 0. 

To first order in exy^ the solution for kx at a magnetic 
interface is 



^ I ^m^r I 1 dz Vv^^ ^m)'^ ^r^rn^xy \ ^ 

V i^m + er) \^ (er " ^my^r^m{^m ^ ^r) ) ^ 

(5) 

Correspondingly, for the decay constants in the y direc- 
tion one obtains the solutions 



l^yd 



er f^u 



i^xyf^x 



%l^yd 



m^xy^x 



(6) 
(7) 



Note in particular the structure of equation (5): the 
right-moving {^{kx) > 0) and left-moving {^{kx) < 0) 
solutions have different magnitudes, because the left- 
right symmetry is broken by the introduction of the 
magnetic field. The symmetry is only restored if the 
sign of the static magnetic field is also reversed. For 
^i^xy) > we have 5R(/c+) > 5R(/c_), and vice versa. This 
'non-reciprocal' behaviour has been observed experimen- 
tally in surface plasmon propagation with a magnetic 
medium^^'^^. Note also how the presence of a real (dissi- 
pative) part in exy^ introducing a non-Hermitian part to 
(1) contributes to the imaginary part of kx and hence in- 
creases the attenuation of the wave, while a purely imag- 
inary exy leaves (1) Hermitian and introduces no addi- 
tional dissipation. 



B. Metal-Insulator-Metal waveguide 

Consider now the geometry shown in figure 1(a). We 
seek to find plasmon solutions, i.e. surface-bound waves 
that in the metal are of the form 

Hm = i^oe^^-^-i^^+'^y-^/2e-'^y-l2^l, > T/2 (8) 

and in the dielectric are linear combinations of decaying 
exponentials from either side 

Hd = {Ae^y^y + Be-^y^y)e'^-''-'^\ \y\ < T/2. (9) 

In general 1^41 7^ \B\ when the dielectric is anisotropic 
(i.e. when exy 7^ 0). 

From the Maxwell equations the electric fields in the 
metal have the form 



_ l^j/m 1^1 rr 
-'-^xm — J^rr 



ujem y 



kx = +^+ or — k- 



(4) 



(10) 



while the £^-fields in the dielectric for an iJ-field of the 
form (9) are 
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^ikxx-iut [ie^k^^-Be-'^y^y - Ae'^y^y) + e^^yi^ydiAe'^y^y - Be-'^y^y)] 

= ^WfTQ ■ 

From the boundary conditions at the two interfaces at = ±T/2 we obtain two equations for A and B. In order for 
them to have non-zero solutions the determinant formed from the coefficients should equal zero. This condition gives 
the dispersion relation for the system: 

{2erem{el + ely)i^ydKym cosh[/^^rfT/2] + (6^(6^^feg + e^^l^) + (e^ + elyfi^U smh[KydT] 
The solutions are 



r^yd cosh[A^^rfT] ± W {^^i^ld + ^lyf^lO^ - cosh[2/^^dT])) /2 



l^yml,2 — T-TT ^ (13) 

In the limit of small Cxy and small cavity thickness T the solution can be approximated by 

-emi^r - ^ly)i^yd tanh[/^^^T/2] em^lyi^l sinh[KydT] 



f^yml 



IC^Kjyd 



-^m{4 - 4y)'^yd coth[A>:^rfT/2] em^^^ smh[/^^^T] 



In the non-magnetic case e 



we recover 
the well-known plasmon dispersion relations K,ymi = 
emtiyd tanh[KydT/2] QQ^Tespondrng to the symmetric {A = 
B) and A>:2ym2 = ^"^^^^ coth[/^^aT/2] ^j^^ antisymmetric 
{A = —5) i^-fields respectively. Note that in the non- 
magnetic limit the Ex field has the opposite symmetry 
to Hz'. Ex is odd if A = B and even if A = —B. In 
this paper we shall always refer to the symmetry (or ap- 
proximate symmetry) of the i^-field when describing a 
mode. 

For large cavities the contribution of cos]i[lkydT] in the 
square root of the solutions of Eq. (13) for the wavevector 
becomes predominant and we recover the single-interface 
solution given by Eq. (7). However now there are 
two forward-propagating and two backward-propagating 
modes: the right-moving mode with kx = /c+ and the left- 
moving mode with kx = —k- are localised on the upper 
interface (which has the same orientation as that consid- 
ered in §11 A), while the other two modes with kx = k- 
and k localised on the lower interface, which 

has the opposite orientation. The additional degener- 
acy between positive and negative k appears because the 
structure is now invariant under a reflection in the line 
^ = 0. 



C. Example system 

To illustrate these effects, we describe a specific exam- 
ple where the magnetic dielectric is bismuth-substituted 



yttrium-iron garnet (Bi:YIG) and the metal is silver. The 
magneto optical data for the Bi:YIG is taken from exper- 
imental data^^'^^. The permittivity for the silver is fitted 
from the Drude model 



= 1 



cj(l + a;2r2) 



(15) 



where ujp is the bulk plasma frequency, with parame- 
ters taken from references 37 and 38. Both permittivities 
are therefore frequency dependent, as shown in Figure 2. 
Throughout we take h = 1 and so give values for angu- 
lar frequency in energy units. The geometry is shown in 
Figure 1(a), which also shows the Ex field distribution 
for a symmetric and anti-symmetric mode. Figure 1(b) 
shows the dispersion relation of the waveguide for differ- 
ent thicknesses when there is no magnetisation in the di- 
electric. The thick blue line is the dispersion relation (Eq. 
3) of surface plasmons on a single interface; for the MIM 
waveguide the mode frequency is split into two modes 
moving in each direction, as expected. The lower lying 
cavity plasmonic modes have a symmetric magnetic field 
profile, while the high energy ones shown with the dark 
red line on the figure are antisymmetric. For thicker cavi- 
ties the energy difference between the symmetric and an- 
tisymmetric solutions gets smaller as one would expect as 
in the limit of infinite waveguide thickness both solutions 
will become almost degenerate, corresponding to the sep- 
arate solutions (Eq. 5) on the upper and lower interfaces 
of the slab. The thick black line is the vacuum light line 
CjO = ckx] solutions to the left of it are inaccessible to 
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a) b) 




1x10^ 2x10^ 3x10^ 4x10^ 5x10^ 
Re[kx] l/[m] 



FIG. 1: (a) The MIM waveguide geometry. The 
coordinate system, showing the direction of the static 
magnetisation and plots of the electric field distribution 
for the antisymmetric (above) and symmetric (below) 
modes — note the symmetry descriptions refer to the 
symmetry of the i^-field (not shown); (b) the solutions 
of the dispersion relation Eq. (12) for as a function of 
kx for the non-magnetic case for a range of cavity sizes. 
The thick blue line is the solution for a single interface. 

The coupled cavity modes are split into low-energy 
symmetric (SM) and high-energy anti-symmetric (AM) 
modes. The black lines are the light lines in free space 
(thick) and in the dielectric (dashed). 



coupling from the outside world in a translationally in- 
variant system. The dashed black line indicates the light 
line in the dielectric medium uj = ckx/s/cd; modes to the 
left of the dielectric light line have propagating solutions 
in the dielectric (i.e. nearly imaginary A^^yd, provided dis- 
sipation in the dielectric is weak). Modes to the right 
of the dielectric light line have nearly real Kyd'', these are 
the bound surface-plasmon solutions that we seek. The 
lowest symmetric mode always has a frequency below the 
plasmon of the semi-infinite surface and remains bound 
for arbitrarily small cavities; the lowest antisymmetric so- 
lution always spills out and becomes a propagating state 
in the dielectric region (i.e. crosses the light line of the 
dielectric) for sufficiently small kx- All modes are rapidly 
decaying into the metallic boundaries of the cavity. 

When a magnetic field is added in the dielectric, 
the dispersion relation is modified. We choose a field 
large enough to saturate the magnetisation, and hence 
the magneto-optic response; in practice 2kOe would be 
sufficients^. For simplicity we give the magnetisation in 
units of the saturation value, so magnetisation along -\-z 
and —z correspond to M = ±1 respectively. The changes 
in kx are shown as a function of energy in Figure 3(a) for 
three different cavity sizes. Modest changes (of the order 
of 1%) are predicted, the sensitivity to magnetisation be- 
ing greatest at high energies where the magneto-optical 
response of the dielectric (figure 3 (b)) is larger. The dif- 
ference in the change of the wavevector for both modes 
is small in the order of 0(e^ ). 




[eV] CO [eV] 




ca [eV] ca [eV] 

FIG. 2: Frequency dependence of the dielectric 
functions: real parts of (a) (b) Cxy] imaginary parts 
of (c) Cr and (d) Cxy. 



Figure 3 (b) shows the wavevector of the propagat- 
ing modes in the x direction for cu = 2.37 eV, which is 
the highest frequency shown on figure 1(b), for waveg- 
uides with different thickness. The wave vectors for the 
single interface are shown for reference. The wave vec- 
tor for the cavity modes doesn't depend on the sign of 
the magnetisation as is clear from the dispersion relation 
Eq.l2. For thin waveguides the symmetric and antisym- 
metric modes are far apart. As the thickness increases 
the wavevector approaches the solutions for a single inter- 
face. The change in the wavevector versus the waveguide 
thickness for change in the magnetisation from M = 
to M = —1 for different energies is shown on figure 3 (c) 
At small thickness the difference in the magnetisation 
modes is insignificant. The waveguides are less sensitive 
for the magnetisation at lower energies. This is due to 
the frequency dependence of Cxy which increases at higher 
energies. The change of the wavevector with the mag- 
netisation is shown on figure 3 (d). For single interface 
the change is linear as is clear from Eq. 5. For cavities 
the dependence is quadratic. The two modes on the two 
interfaces couple to produce the anti-crossing. The mag- 
netic field changes further the levels. For small thickness 
the change of the wavevector due to the coupling of the 
two interfaces is predominant. At larger waveguide thick- 
ness the coupling between the interfaces weakens which 
makes the magnetic field effect more pronounced. In the 
limit of large cavities the solution will follow the one for 
a single interface. 

The intensity of a surface plasmon propagating along 
an interface decays as ex.p{—2^[kx]x) , owing to the en- 
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FIG. 3: (a) The change in wavevector Sk = k-^ — k 
resulting from magnetisation of the dielectric as a 
function of frequency for different cavity thicknesses; 

(b) dependence of wavevector on the thickness of the 
waveguide for uj = 2.37 eV and different magnetizations; 

(c) absolute change in the wavevector vs. waveguide 
thickness when the magnetization changes from M=0 to 

M=-l for different energies ; (d) dependence of the 
wavevector on the magnetisation for different cavity 
thicknesses at = 2.37 eV; 



ergy loss from Joule heating; the corresponding prop- 
agation length is Lp = 2^[fc ' Figure 4 (a) shows 
the propagation length for symmetric and antisymmet- 
ric modes for different cavity thickness. Generally the 
propagation length rises as the frequency is reduced, es- 
pecially for modes which approach or cross the dielectric 
light line so the fields become less strongly bound to the 
metal-dielectric interfaces; the plasmons can then propa- 
gate long distances of more than 20 /i m. However, once 
solutions reach the minimum frequency for propagation 
along the waveguide without penetrating the electrodes 
(i.e. approach the frequency axis in Fig.l), the fields 
are forced to penetrate the metal once again as the fre- 
quency is reduced further and the propagation lengths 
drop rapidly. For a given type of mode the propagation 
length Lp becomes shorter as the cavity thickness T is 
reduced. 

In figure 4(b) we show the normalised propagation 
length given by Lp/Xp which is important if we want 
to apply these waveguides as cavities as it gives the num- 
ber of oscillation cycles before decay and is therefore a 
key factor in the determination of the sharpness of the 
resonance i.e. the quality factor of the cavity. The anti- 
symmetric modes for small cavities have much less than 
one cycle of oscillation for smaller energies where there 
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(a) (b) 




1.2 1.4 1.6 1.8 2.0 2.2 1.2 1.4 1.6 1.8 2.0 2.2 
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FIG. 4: Decay lengths for a range of cavity sizes 
(non- magnetic case) as a function of uj: (a) the 
propagation length Lp; (b) the normalized propagation 
length ^{kx)Lp/ {27r); (c) the skin depth in the 
dielectric S^] (d) the skin depth in the metal Sm- 



are no propagating solutions in the waveguide, in accor- 
dance with the dispersion relation. For thicker cavities 
and symmetric modes there are several cycles of oscilla- 
tions, but the oscillations decay rapidly for energies above 
2 eV. 

Figures 4 c) and d) show the skin depths in the di- 
electric core in the metal claddings and Sm- The skin 
depths are defined as 6d = l/^[tiyd] and Sm = l/^[tZym] 
respectively. Note that the the smaller skin depth at 
larger energies of the incident light doesn't correspond 
to longer propagation length; the losses are instead de- 
termined by the proportion of the field energy forced to 
reside in the metallic regions, which is in turn determined 
by the geometry of the cavity. 

Fig. 5 shows the propagation length for the symmet- 
ric and anti-symmetric modes for a waveguide with a 
thickness of 200 nm and for semi-infinite interfaces for 
different magnetizations. For the symmetric mode the 
magnetization increases the propagation length while for 
the anti-symmetric mode it reduces it. This behaviour is 
implied by Eq. (13), which indicates that the magnetisa- 
tion has different effects on the wavevectors of different 
modes. The difference between the propagation length 
for different magnetizations for the waveguide modes is 
insignificant for low energies but becomes larger than the 
corresponding difference for a single interface at higher 
energies. For a semi-infinite interface the propagation 
length is different for the positive and negative directions 
of the magnetisation, as expected from the dispersion re- 
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lation given by Eq. (5). 



III. MAGNETIC CONTROL OF THE LIGHT 
FIELD 

In a non-magnetic plasmonic waveguide the fields have 
a maximum value at the interfaces and are exponentially 
decaying away from them. The field has a symmetric or 
antisymmetric profile inside the waveguide, and the en- 
ergy density is symmetrically distributed along the two 
interfaces. In the case of a magnetic waveguide the mag- 
netic field will perturb the symmetry. To illustrate the 
effect of these differences, and to show how the relatively 
small magnetic modifications to the dispersion relation 
can nevertheless have large effects on the field distribu- 
tions, we employ finite-element method simulation us- 
ing commercial software (COMSOL Multiphysics) to find 
the field distribution in a BiiYIG magnetic waveguide of 
80 nm thickness and length 400 nm confined between two 
silver plates. We excite the radiation by a point dipole 
emitter placed in the centre of the structure. We work 
at two different excitation frequencies: at the higher fre- 
quency uo = 2.37 eV the waveguide has two propagat- 
ing plasmonic modes (one approximately symmetric and 
one approximately antisymmetric) while at the lower fre- 
quency uj = 2eV there is only one such mode (approxi- 
mately symmetric). 

Figure 6 shows the electric and magnetic fields and 
the energy distribution as a function of y (i.e. distance 
from the centre line of the dielectric) for the MIM struc- 
ture, 10 nm from the end, for excitation at the higher 
frequency uj = 2.37 eV. At zero magnetisation M = 
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FIG. 5: The propagation length for cavity of thickness 
200 nm for different magnetizations. 



the fields are completely anti- symmetric (for H^) or sym- 
metric (for Ex). This refiects the dominant excitation of 
the mode even in the electric field by the centrally placed 
dipole emitter. The magnetisation breaks the symmetry 
and as a result the zero of the Hz field is shifted from 
the middle of the waveguide. The energy density is still 
localised along the interfaces but in the presence of mag- 
netisation it is localised preferentially along one of them. 
The maximum field amplitude is correspondingly differ- 
ent at the two interfaces. 

Figure 7 shows the corresponding field patterns for ex- 
citation at the lower frequency {uj = 2.0 eV). In the case 
where M = symmetry forbids any coupling to the mode 
antisymmetric in , so all the radiation is absorbed into 
the mode even in Ex (odd in i^^); because this mode can- 
not propagate, the fields near the ends of the cavity are 
very small. In the presence of magnetic field the symme- 
try is broken and the zero of the Ex-fi.e\d of the propagat- 
ing mode is shifted from the centre, allowing coupling to 
the emitter and strong radiation into the mode. Hence 
for M 7^ the dominant excitation is of the mode nearly 
symmetric in Hz (nearly antisymmetric in Ex): although 
this mode cannot couple to the emitting dipole at M = 0, 
it is the only one that can propagate along the waveguide 
at this frequency. 

This shift in the node of the electric field across the 
cavity with the magnetisation allows us to control the 
coupling of the dipole to the cavity modes and effectively 
switch it on and off. To demonstrate this we simulate an 
oscillating electric dipole of magnitude d with a current 
dipole moment ujd = 1 nA and frequency uj = 2.37 eV, 
oriented along the x axis and placed in the middle (at 
coordinates x = and = 0) of a cavity which is 400nm 
long and 40 nm thick. Such a narrower cavity supports 
only one propagating plasmonic mode (symmetric in Hz 
and anti-symmetric in Ex) even at this higher frequency. 

Figure 8 shows the Ex field for different magnetiza- 
tions. At zero magnetisation (M = 0) the electric field 
is exactly zero in the middle of the waveguide for the 
propagating mode, hence this mode is not excited. The 
electric dipole excites only the antisymmetric (i.e. sym- 
metric in Ex) mode which is very rapidly decaying; the 
resulting field is strongly bound to the central region of 
the waveguide. When M 7^ 0, as was shown on figure 7, 
the electric field is also different from zero at the posi- 
tion of the dipole for the symmetric (i.e. antisymmetric in 
Ex) mode; hence both modes are excited. The symmetric 
mode has a significantly longer propagation length and 
even though little energy couples into it, it is carried away 
across the interfaces. The antisymmetric mode decays 
very fast and at less than half of a plasmon wavelength 
away from the dipole its electric field becomes equal to 
that of the symmetric mode. This leads to a relatively 
complex interference of the two modes in this region. 

These effects also have consequences for the radiation 
which couples out of the cavity. We study the radiation 
emitted into free space for the case where the cavity is 
40 nm thick but only 200 nm long and where the silver 
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FIG. 6: Field and energy distributions as a function of distance from the waveguide median for a Ag/Bi:YIG/Ag 
structure with core thickness T = 80 nm at frequency uo = 2.37 eV: (a) the electric field E^'^ (b) the magnetic field 

Hz] (c) the average total energy density. 
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FIG. 7: Field and energy distributions as a function of distance from the waveguide median for a Ag/Bi:YIG/Ag 
structure with core thickness T = 80 nm at frequency uj = 2 eV: (a) the electric field ; (b) the magnetic field ; 

(c) the average total energy density. 
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FIG. 8: The Ex field from a radiating electric dipole 
within a cavity of thickness 40 nm and length 400 nm 
for different magnetisations of the dielectric 



of emission is dominated by the geometry of the slot and 
its effect on the relative phase of the fields around it. 

In Figure 10 we show the emitted radiation in the far 
field, showing several different geometries in order to il- 
lustrate the different regimes. The first three subfigures 
are for a 400nm long cavity. At = 2.37 eV (figures 10 a) 
and b) the radiation at zero magnetisation is largely sup- 
pressed, since (as shown previously) no energy couples to 
the symmetric propagating mode (anti-symmetric m. E^)^ 
while the antisymmetric mode has very short propagation 
length. In the presence of magnetisation the electric field 




layers are 300 nm thick. The resulting £^a^-field is shown 
on figure 9 for different magnetizations at distances up to 
1 fiia around the structure. As can be seen, the magnetic 
field skews the emitted electric field pattern: the angle 



FIG. 9: The near-field distribution of Ex outside a 
cavity of thickness 80 nm thick and length 200 nm 
containing a radiating dipole at the centre, for different 
magnetisations of the dielectric 
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a) b) c) 




FIG. 10: Polar plots of the far-field radiation norm from a dipole placed in the centre of a cavity filled with 
magnetic or non-magnetic BiiYIG. The magnetic state of the dielectric is given by the colour or the curves: (a) for a 
cavity of thickness 40 nm and width 400 nm at = 2.37 eV (b) for a cavity of thickness 80 nm and width 400 nm at 
uo = 2.37 eV (c) for a cavity of thickness 80 nm and width 400 nm uo = 2eV. (d) for a cavity of thickness 40 nm and 
width 200 nm at = 2.37 eV (e) for a cavity of thickness 80 nm and width 200 nm at = 2.37 eV (f) for a cavity of 

thickness 80 nm and width 200 nm at = 2 eV 



distribution is only slightly different for different mag- 
netisation; the main role of M is to break the symmetry 
and turn on the radiation. For shorter cavities, shown in 
figures 10 d) and e), the radiation emitted into free space 
at zero magnetisation is due to the energy coupled to the 
antisymmetric mode. When the magnetisation is differ- 
ent from zero the two modes are both excited; although 
they have different degrees of excitation in the centre of 
the structure, they decay so that the field amplitudes at 
the edge are comparable. The magnetic field then tunes 
the relative phases of these contributions and alters the 
interference from constructive to destructive, depending 
on the direction. The resulting interference pattern of the 
far- field radiation is therefore skewed to different sides 
depending on the magnetisation. Finally for 80nm cav- 
ity with a dipole oscillating at = 2 eV in short cavities 
such as the 400 nm (figure 10 c) ) and the 200 nm (figure 
10 e) ) the antisymmetric mode has not decayed by the 
end of the cavity; hence the radiation patterns at M = 
and for M 7^ are similar. 



IV. CONCLUSION 

The modes of a magnetised MIM waveguide form an 
interesting contrast to the modes of a single magnetic in- 
terface. Because the system as a whole has a plane of 



symmetry at = 0, the allowed values of the wavenum- 
ber are not direction dependent for a given frequency. In- 
stead, the modes come in pairs with very different weights 
of the electric field on the upper and lower waveguide sur- 
faces. 



Even though the changes in the dielectric function in- 
troduced by the presence of the magnetisation may be 
small in absolute terms, these qualitative changes in 
the nature of the optical modes can lead to significant 
changes in the response of the structure to excitation. 
We have illustrated this by computing the field distri- 
butions in a finite-length slot waveguide in response to 
excitation by a point dipole placed at its centre. Both 
the total emission from the cavity and the pattern of far- 
field radiation can be strongly modified by switching the 
magnetic field; the total emission is predominantly con- 
trolled by the presence or absence of propagating modes 
with strong coupling to the radiating dipole, while the 
far-field radiation pattern is controlled by the phase dif- 
ferences at the opposite ends of the guide. These exam- 
ples point to the possibility of using magnetic control to 
switch the propagation of fields in more complex photon- 
ics structures. 
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